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Chapter 1 
Introduction 


1.1 Motivation 

Le hasard est un phenomene qui existe dans la nature, alors comme tous 
les autres phenomenes naturels, on veut le comprendre et l’etudier d’une 
maniere scientifique, le seul probleme c’est que le hasard, par definition, est 
non sounds a aucune une loi, par exemple, lorsqu’on lance une piece de mon- 
naie, on ne peut pas avoir des lois scientifiques qui nous aide a determiner 
d’une maniere certaine le resultat qu’on va obtenir avant de faire le lance- 
ment. done pour mieux comprendre le hasard, l’une des methodes porpses 
est de le simuler, alors ca veut dire quoi une simulation? 

Revenons a l’experience de la piece de monnaie, une simlation de cette 
experience sert a trouver une methode pour obtenir des resultat tres sim- 
ilaire aux resultats qu’on obtient en realite, mais sans etre oblige de faire 
1’experience en realite, ceci est possible grace aux etude des probability qui 
nous donnent des infomations importante sur les phenomene aleatoires qui 
vont nous permettre de les simuler, par exemple, au point de vue proba- 
biliste, on salt que l’experience de la piece de monnaie suit une loi uniforme 
de parametre 1/2, alors le but devient simuler un variable aleatoire uniforme 
de parametre 1/2. Alors le but finale de ce document est de developper un 
programme qui fait la simulation des variables aleatoires reellc usuelles. 

1.1.1 Pourquoi faire des simulations? 

La simulation des variables aleatoires est utilisee dans beacoup d’applications, 
par exemple: 
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• Etude des phenomenes aleatoires: la simulation d’une loi d’une vari¬ 
ables aleatoires nous permet d’obtenir un echantillion statistiques qui 
nous permet de comprendre tous les phenomenes aleatoires aui suit 
cette loi, mais pour etudier ce phenomene on a besoin de faire des mille 
ou meme des millions de simulations, ce qui est tres couteux ou n’est 
meme pas possible dans la vie reelle. 

• Etude des sciences naturels: les phenomenes physiques, geologiques, 
biologiques, chimiques, sont tous influence par le hasard,alors pour bien 
comprendres ce qui passe dans ces phenomenes, les savants ont besoin 
de faire des simulations de variables aleatoires. 

• Developpment des jeux et des envirenments virtuclle: la simulation des 
variables aleatoires permet de developper des jeux et des envirenment 
ou les phenomenes aleatoires comportent comme ceux de la vie reelle. 


1.2 Plan du document 

Ce document est divise en deux parties, une partie therorique est une 
partie pratique, dans la partie theorique, on va parler des variables 
aleatoires qu’on veut simuler,puis on va demontrer la methode de simu¬ 
lation de chaqune de ces lois, et puis on va parler des tests d’hypotheses 
comments les utiliser pour verifier nos simulations, et puis, dans la par- 
tie pratique, on va utiliser utiliser les resultats qu’on a etabli dans la 
partie pratique pour developper un programme qui fait la simulation de 
chaqune des loi qu’on a vu, et puis on va appliquer les test d’hypotheses 
pour verifier que la simulation de notre programme est bonne. 
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Part I 

Partie theorique 


6 



Chapter 2 


Rappels: Elements de 
probabilites 

2.1 Variables aleatoires discretes 

2.1.1 Experience aleatoire 

Definition 2.1.1. On appelera une experience aleatoire toute experience 
dont on ne peut pas prevoir lc resultat de fagon certaine, void quelques 
exemples: 

Exemple 2.1.1. On lance un piece de monnaie sur une table cinq fois, et 
on prend comme resultat le nombre de fois ou on obtiendra une ’face’. 

Exemple 2.1.2. On lance un de et on note comme resultat le nombre obtenu 
sur la face superieure du de. 

2.1.2 Variable aleatoire 

Exemple 2.1.3. Considerons le jeu du lance d’un de, les resultats possibles 
sont: 1, 2, 3, 4, 5, 6. Notons 0 V ensemble de tous les resultats possibles, 

Q = {1,2,3,4,5,6} 

chaque fois qu’on lance le des, on obtient comme resulat un element de 0, 
supposons que le resultat obtenu est 1, on dira que l’evenement : ’obtenir 
1 comme resultat’ a ete realise, on peut s’interesser a des evenements plus 
complexes qu’un simple resultat, par exemple: 
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A = ” Obtenir un nombre paire” 

B = “lc resultat est un nombre plus grand qne 4” 

On dira que l’evenement A est realise chaqne fois qu’on obtient un nombre 
paire. ces evenements penvent etre ecrit d’une maniere explicite: 

A = {2,4, 6} i? = {5,6} 

et sont tous des sous ensembles de 0. 

Notons par to lc resultat de L experience aleatoire , alors d’une maniere gen¬ 
erate, on dira que A s’est realise si u E A. 

Definition 2.1.2. On appelera l’espace fondamental ou univers, l’ensemble 
de tout les resultats possibles de l’experience aleatoire. et on le not par 0 . 

Definition 2.1.3. On appelera tout sous ensemble de O un evenement. 

Definition 2.1.4. On appelera une probabilite toute application P definit 
par: 


Et qui verifie: 


rm-M o,i] 

\ A P(A ) 


• P(Q) = 1 

. Vi, j e {1, 2, ...n}/Ai O Aj = 0 : ULi p ( A k) = ELi p ( A k) 

Definition 2.1.5. On appelera une variable aleatoire reelle tout application 
X definit par: 


jn X(Q) C R 
|a4X(A) 

Exemple 2.1.4. On lance un de deux fois et on pose: 


X = {La somme des deux resultats obtenus}. 



Notation 


Soit X une variable aleatoire, u un evenement de 0, et a € M. 

Dans le rest de ce document, on va suivre les notations suivantes: 

P(X ^ a) := P({u £ fl/X(u) ^ a}) = - oo, a]) 

P(X = a) := P({co £ Q/X(cj) = a} = ^({a}) 

Definition 2.1.6. Soit X une vairbale aleatoire reelle, on definit la fonction 
de repartition de X, notee Fx, par: 

fA'(fi) -> [0,1] 

\ P(X < x) 

Definition 2.1.7. Une variable aleatoire X est elite discrete lorsque l’ensemble 
des valeurs qu’elle peut prendre est denombrablc. 

Exemple 2.1.5. Dans l’exemple precedent, la variable aleatoire X definie 
par: 


X = {La somme des deux resultats obtenus}. 
est une une variable aleatoire discrete. 

Definition 2.1.8. Soit X une variable aleatoire reelle discrete, l’esperance 
de X, notee par E(X) est definie par: 

E(X) = Etexm kP(X = k) 

Definition 2.1.9. Soit X une variable aleatoire reelle discrete, la variance 
de X, notee par V(X) est definie par: 

V(X) = E(X 2 ) - E(X) 2 

Interpratation de l’esperance et la variance 

En probabilites et statistiques, l’esperance et la variance nous donnent des 
informations importantes sur un echantillion statisticpies ou un phenomene 
aleatoire. L’esperance represente la moynenne des valeurs d’une variable 
aleatoire A", e’est pour ceci qu’on l’appelle aussi la moyenne, prenons un 
exemple pour bien comprendre: supposons par exemple que X prend des 
valeurs dans l’ensebmle des points des etudiants en examen note sur 20, 
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apres les calculs, on trouve par exemple que l’esperance est egale a 14 et que 
la variance egale a 2, ceci signifie que la note moyenne des etudiant dans 
cet examen est 14,et que la variance des notes ne depasse pas 2, c’est a dire 
que les notes de la plupart des etudiants dans cet examen varient entre 12 
et 16, done l’eserance et la variance nous donnent une tres bonne idee sur 1c 
comportement de la variable aleatiore X. 

2.1.3 Loi de probability 

Definition 2.1.10. En theorie de probabilities, une loi de probability decrit 1c 
comportement d’une variable aleatoire X, elle est definie souvent en utilisant 
Fesperance et la variance de X. 

2.1.4 Loi de Bernoulli 

Soit un univers const it ue de deux evenements: 

n = {v,F} 

Sur 11 on definit une variable aleatoire X par: 

X(V) = p X{F ) = 1 -p p^ 0 
alors : X(H) = {0,1} 

Definition 2.1.11. On dit que X suit une loi de Bernoulli de parametre p, 
et on note : 


X ~ B(p) 


2.1.5 Loi binomiale 

Soit Xi, X 2 ,X n n variables aleatoires suivant une loi de Bernoulli de 
meme parametre p6 [0,1], soit X une variables aleatoires definit par: 

X = X 1 + x 2 +.... + x n 

alors: X(H) = {l,...,n} 

Definition 2.1.12. On dit que X suit une loi binomiale de parametres n et 
p, et on note : 

X ~ B(n,p) 
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2.1.6 Loi de Poisson 

Soit X une varibale aleatoire qui verifie: A"(f2) = N 

Definition 2.1.13. On clira que X suit une loi de Poisson de parametre A 
si: 

VA; e N P(X — k) — 

2.2 Variables aleatoires continues 

Soit X une variable aleatoire reelle. 

Definition 2.2.1. On dira que X une variable aleatoire reelle continue,ou 
bien a densite, si il existe une fonction reelle fx, qui verifie: 

Va, b E M, a ^ b P{a ^ X ^ b) = J c( 5 fx(x)dx 

fx est appelee la densite de X. 

Remarque. Remarquons que: 

Va e M F x (a) = fx(x)dx 

Definition 2.2.2. On definit 1’esperance d’une variable aleatoire reelle con¬ 
tinue par: 

E ( x ) = IZo x fx{x) dx 

2.2.1 Loi uniforme 

soit a,b G M 

Definition 2.2.3. On dit que X suit une loi uniforme sur le domaine [a,b\, 
si X est une variable a densite fx definie par: 

—y M 

fx : < x ^ ^ si x e [a, b] 

' x (->• 0 sinon 

Et on note : 

X~U([a,b]) 


11 



2.2.2 Loi exponentielle 

Definition 2.2.4. On clit que X suit une loi exponentielle sur le domaine 
de parametre A, si X est une variable a densite fx definie par: 

/■M —y M 

fx - < x i-» \e~ Xx si x ^ 0 
^ x i —y 0 sinon 


Et on note : 


X ~ £(X) 


2.2.3 Loi normale 

Definition 2.2.5. On dit que X suit une loi normale sur de parametres m 
et a, si X est une variable a densite fx definie par: 

(M -> M 

IX- \ 1 -1 / x-m -.2 

1 -7==- e 2 i > 

l vz-rrcr 


Et on note : 


X ~ A f(m, a) 

2.2.4 Lois aux densites de support compact 

On dit que X suit une loi normale sur de parametres m et a, si X est une 
variable a densite fx de support compact, i.e: 

3a, b e M : Supp(f) = {ieR / fx(x) 0} = [a, b] 
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Chapter 3 

Simulation des variables 
aleatoires : Partie theorique 


Dans la plupart des languages de programmation, on trouve une fonction 
predefini qui retourne apres chaque appcl une valeur reclle aleatoire entre 0 
et 1 (par exemple rand dans le language C), au point de vue probabiliste, cett 
fonction suit une loi uniforme sur le domaine [0,1], lc but de ce chapitre et de 
trouver un algorithme qui perrnet d’utiliser cette fonction pour simuler toutes 
les lois qu’on a vu, au point de vu mathematique, ceci revient a etablir des 
relations mathematiques entre la loi uniforme sur lc domaine [0,1] et toutes 
les autres lois usuelles. 

Pour generaliser les notions, on appelcra cette fonction un generateur, note 
par U, et a chaque fois qu’on appellc cette fonction, on dira que ’On tire 
le generateur U\ cette expression signifie au poit de vu de programmation 
qu’on fait appel a une fonction qui retourne une valeur reellc aleatoire entre 
0 et 1. 

Dans le reste de ce chapitre on note par X la variable aleatoire qu’on veut 
simuler. 


3.1 Loi uniforme 

On veut simuler une loi uniforme dans un domaine [a, b], puisqu’on a deja 
un generateur suivant une loi uniforme qui retourne un valeur u entre 0 et 1, 
il suffit alors de centrer cette valeur dans le domaine [a, b], done l’algorithme 
suivi est: 
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• On tire U, et on obtient une valeur u 

• On prend x = a + (b — a)u 

Par conseqence la valeur de x est une valeur aleatoire sur le domaine [a, b], 
et on ecrit X = a + (b — a)U. 

Proposition 3.1.1. La variable aleatoire obtenue X suit une loi uniforme 
sur [a,b]. 

Preuve. Soit c G [a, b], on a: 


P{X ^ c) = P(a + (b- a)U < c) 


= P{U < 


c — a 
b — a 


(3.1) 


Comme c G [a, b], alors : 


0 < ^ i 

^ b—a 


Done: 


P(U = s 


Alors: 


F x (c) = P(X^c) = fEf Vc G [a, b] 


Done X suit une loi uniforme sur [a,b]. 


3.2 Loi de Bernoulli 

Supposons que X suit une loi de Bernoulli de parametre p G [0,1], et notons 
que: 

• P(X = 1 )=p 

• P{U ^ p) = p 

Alors les deux evenements {A" = 1} et {U ^ p} ont 1a. meme probability 
done l’evenement (A" = 1) peut etre simule par l’evenement (U <= p), et 
done l’algorithme suivi est: 
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• On tire U, on obtient line valeur u comprise entre 0 et 1. 

• si u <= p, on prend x — 1 , sinon, on prend x — 0 . 

Alors X = t{U ^ p} 

Proposition 3.2.1. La variable aleatoire obtenu X suit une loi de bernoulli 
de parametre p. 

Preuve. On a: 

A"(f2) = = {0,1} 

Et on a: 

(X = 1) = (1 {U ^p} = 1) <*=*► U < p 

Alors: 

P{X = 1) = P(U < p) 

Comme U suit une loi uniforme sur [0,1] et p G [0,1] alors: 

P(U ^ p) = p 

Done: 

P(X = 1 )=p 

Et on a : 

(X = 1) <*=> 1{U < p} = 0 <*=*► (U>p) 

Alors: 

P(X = 1) = P(U >p) = 1 - P(X ^p) = l-p 
Done X suit une loi de Bernoulli de parametre p. 
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3.3 Loi binomiale 


Supposons que X suit une loi binomiale de parametres n et p, X rv_/ B(n; p). 

On sait qu’une variable binomiale represente la somme de n variables independantes 
de Bernoulli de parametre p, alors pour simulier la varibale X, il suffit de 
simuler n variables aleatoires independantes suivant une loi de Bernoulli et 
de prendre la somme des valeurs obtenu comme resultat, alors: 

-X = YJk =1 Y k avec ■ (be) ~ B(p) v/c G {1, n} 

En utilisant le resultat etabli pour la loi de bernoulli, ceci est equivalent a: 

-X = Y2 =i 1 {u k ^p} avec '■ U k ~ U([0,1]) V/c G {1, n} 

Proposition 3.3.1. La variable aleatoire obtenue X suit une loi binomiale 
de parametres n et p. 

Preuve. Comme la variable aleatoire l{n<p} Sll it une loi de Bernoulli de 
parametre p, alors d’apres la definition d’une loi binomiale, la variable X 
definit par: 

X = Y!l= 1 1 {u k ^ P } avec- : U k ~ U([0,1]) VfcG {1, n} 
suit une loi binomiale de parametres n et p. 


3.4 Loi de Poisson 

Supposons que X suit une loi de poisson de parametreA: 


X ~ P(A) ^ Pk = P(X = k) = e~ A * ^ 
On note par P k le cummul des p k : 

Pk = Ej=i Pk 

Proposition 3.4.1. La variable aleatoire X definie par: 

X = £* i ^U< P k } 

suit une loi de poisson de parametre A. 
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Preuve. Supposons que : 


Montrons que: 


X = En=1 HPk-l <U< P k } 


X ~ P(A) 


On a: 


P(X = k) = P(P fc _l < U<P k ) 


Et notons que: 

V/c G N Pfe+i = Pfc + 

Alors: 

P(X = k) = P(P fc _! < U < P k —\ + ^ fc _0 

Done: 


P(x = k) = P(0 ^ U < ^p k -i) 

Comme U est une loi uniforme sur [0,1], alors: P(0 ^ U < u) = u Done: 

P(X = k) = = e- A f 

Alors: X suit une loi de Poisson de parametre p. 

Alors si X est une loi de poisson de parametre p, on peut la simuler en 
ntilisant l’expression: 


X = En=1 H p k -1 <U< P k } 

Remarque. Au niveaux de programmation, cette methode n’est pas pratique 
car on ne peut pas simuler l’infini pratiquemet, la solution qui reste est de 
simuler l’infini par une tres grande valenr, plus que la valeur est grande, 
plus que le resultat obtenu est plus precis, rnais lorsqu’on choisit une valeurt 
grande, la machine doit effectuer trop de calculations avant de determiner 
le resultat, ce qui cause le programme d’etre lent, heureusement pour la loi 
de Poisson, il y’a une autre methode qu’on va voir qui resoudra tous ces 
problemes. 
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3.5 Loi exponentielle 


Pour simuler la loi exponontielle, on a besoin de theoreme suivant: 

Theoreme 3.5.1. Soit X une variable aleatoire de fonction de repartition 
F strictement croissante, alors on a: 

F(X) suit une loi uniforme sure le domaine [0,1], i.e : F(X) r\j U([0,1]) 

Preuve. Soit X une variable aleatoire de fonction de repartition F strictement 
croissante, Par definition on a : 

F(x) = P(X ^x)e [0,1] 


On pose : 


u = F(x) 


x = F 


-i/ 


u) 


On remplace: 


F{F~\u)) = P{X < F-\u)) 


Done: 


u = P(X ^ F~\u )) 

Comme F est strictement croissante, alors: 

X < F~\u)) <*=*► F{X) < u 

Done: 


u = P(F(X) ^ u ) 

Ce qui montre que: F(X) ~ U([0,1]) 

Si la variable aleatoire X admet un fonction de repartition strictement 
croissante, alors on a montre que : F(X) ~ U([0,1]), ce qui est equivalent a 
dire que: F _1 (C7) suit la loi de X. 

Done si on commit la fontion inverse F -1 , pour simuler la variable X, il suffit 
alors de prendre X = F~ 1 {U ), la loi de Poisson est une application de cette 
methode: 

Exemple 3.5.1. X ~ £(\) F(x) = 1 — e~ Xx =^- F\x) = \e~ Xx 

Alors : Vr G R F'(x) > 0 F est strictement croissante. 

Done d’apres lc resultat etabli, pour simuler la loi expnontiellc, il suffit de 
prendre: X = F~\U ) <*=*► X = 

Mais puisque U suit une loi uniforme sur [0,1], alors 1-U suit la meme loi, 
done on peut prendre: X = — 1 "] 1 ' J . 
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Relation entre la loi de poisson et la loi exponentielle 

Si des si des evenements surviennent a des dates separees, alors la duree entre 
deux evenemets suit une loi de exponentielle d’un parametre, et le nombre 
d’evements survenant dans une unite de temp suit une loi de poisson du meme 
parametre. Par exemple la duree entre lc passage de deux voiture dans une 
autoroute suit une loi exponentielle, et lc nombre de voiture qui passe dans 
un unite de temps suit une loi de poisson. Alors si on peut simuler des vari¬ 
ables aleatoire exponentielles de parametre A, on peut obtenir une variable 
aleatoire de poisson du meme parametre en comptant lc nombre de simula¬ 
tions necessaires pour que la somme des valeurs des variables aleatoire expo- 
nentielles simules depasse la valeur 1. II est facile a voir que cette methode 
et plus pratique que la methode precedente au niveaux de programmation. 

3.6 Lois aux densites de support compact 

Supposons que X est une vairable aleatoire de densite / a support compact 
[a, b] avec a,6eK Alors : 

3M > 0 telque : VieR |/(x)| ^ M 



On a: 
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• surface du rectangle (abed) — M(b — a) 

• surface sous la courbe = 1 (puisque e’est une densite) 

• surface de la zone hachuree = f a '° f(x)dx = f(x)dx = F(x 0 ) 

Pour simulcr la varibale X, on va utiliser une methode qu’on appelle ralgorithme 
de rejet: 

• Sur la surface du rectangle abed on choisit d’une maniere aleatoire 
uniforme un point A(X, Y), pour faire ceci, on applique le precede 
suivant: 


• On tire un generatenr U\ et on pose: X — a+ (b — a)Ui 

• On tire un generatenr t/ 2 et on pose: Y = Mf/ 2 


• Si Y ^ f(X) on garde X 


• Sinon, on repete le precede jusqu’a ce qu’on obtient: Y ^ f(X ) 

Theoreme 3.6.1. La varibale aleatoire obtenue ’X sachantY ^ f(X)” suit 
une loi de densite f. 

Preuve. Pour montrer quo la variable ” X/Y ^ /(X)” suit une loi de densite 
/, il suffit de montrer que: 




Sachant que: Y ^ f(X), soit xq G M Alors : 

P(X n Xo/Y ^ /(A')) = 


On pose: 

• Sh = Surface hachuree 

• Sabcs = Surface du rectangle abed 

• Sf = Surface sous la courbe = 1 

p(x^ Xo nY^f(x)) = 

P (Y < /(A)) = ^ 


s h _ F(xp) 
Sabcd. M(b—a) 


1 

M(b—a) 


Alors : 


P(X ^ x 0 /Y ^ f(X )) = F(x 0 ) 
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3.7 Loi normale 


La loi normale n’a pas line densite a support compact et on ne connait pas 
d’expression simple de Finverse de sa fonction de repartition. On ne peut 
done employer aucune des deux methodes precedentes. 


3.7.1 Methode de Box-Muller 


Pour simuler la loi normale on doit appliquer la methode de Box-Muller qui 
perrnet de simuler deux variables aleatoires nor males centrees reduites inde- 
pendantes : 

Soit X et Y deux variables aleatoires normales centrees reduites indepen- 
dantes, i.e : X ~ N(0,1) et Y ~ N(0,1) 

Pusique X et Y sont independantes, alors la densite jointe de X et Y: 


fx,r(x,y) = fx(x)f Y (y) = h e " 2 


On effectue le passage en coordonnees polaires: 


x — r cos(0) y — r sin(6 ) ) 


En appliquant le theoreme de changement de variables, on obtient: 


.I'x,y( x , y)dxdy = fx,Y( rcos (^)i r sm(9))rdrd9 

= rdrdO (3.2) 

27T 

= /i?,0(r, 0)drdQ 


fR,e represente la densite jointe des variables R et O, dans cette densite on 
recommit: 

= densite de la variable O, alors O suit une loi uniforme sur [0, 27 t]. 

-r 2 

re 2 = densite de la vairiable R. 


On en deduit la fonction de repartition de R : 

-t 2 -r 2 

F R (r) = Pr{R ^ r) = f Q r te~dt — 1 — e~ (Integration par parties) 
Alors: 
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F R 2 (r) = Pr(R 2 ^ r 2 ) = Pr(R ^ r) = 1 — e 2 
On pose: u = r 2 et on obtient: 

Fr 2 (w) - 1— r: 2 " 

Alors on trouve que R 2 suit une loi exponontielle de parametre \ : 

R 2 ~£(\) 0~U[O,2vr] 

Done d’apres les resultats deja etablis, si on tire deux generateurs U\ et U 2 , 
la simulation de M et 0 est donnee par: 

R = \J (—2ln(Ui)) 0 = 2nU 2 

Remarque. On salt que si A" ~ A/”(0,1), et si Y — rn + crA" alors: 

Y ~ A f(m, cr) 

Done pour simuler une variable Y ~ A f(m, cr) on sirnule X ~ A/"(0,1), et on 
prend: 

Y — m + aX 

3 . 7.2 Methode de limite centrale 

Une aute maniere qui perrnet de simuler la loi normale est une application 
de theoreme central limite, enonce ci-dessous: 

Theoreme 3.7.1. Soit X 1 , X 2 ,X n n variales aleatoires independantes 
suivat la meme loi d’esperance m et de variance a, et soit S n leur somme : 
S n — Xi + ... + X n alors la variable aleatoire definie par: 

Sn—nm 

y/na 

converge vers une loi normale centree reduite lorsque n tend vers l ’infini. 

Alors cette methode consiste a tirer un grand nornbre n de generateurs 
Ui,U n et puis prendre la somme Z n — U\ + ... + U n , et done la varibale : 



est tres proche d’une loi normale centree reduite. 

Remarque. En effet, la convergence du theoreme central limite est tres vite, 
on peut obtenir des bonnes approximations pour n — 10, mats la methode de 
Box-Muller reste plus pratique car elle donne des simulations exactes, meme 
au point de vue la machine, la methode de limite centrale demande plus de 
caluclations que la methode de Box-Muller. 
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Chapter 4 


Test d’hypotheses 


4.1 Introduction 

Un menuisier a commande a un fabricant 500 plaques de taille de 150cm, 
apres avoir regu sa commande, le menuisier se demande si, a cause d’une 
erreur, lc fabricant lui a envoye des plaques de taille 155cm au lieu de 150cm, 
pour s’assurer qu’il a obtenu la bonne commande, lc menuisier prend un 
echantillion de 50 plaques d’une maniere aleatoire, mesure chaque plaque, 
et calcule la moyenne des mesures, d’une maniere intuitive, lc menuisier 
a decide que si la valeur de la moyenne est entre 149cm est 151cm, alors 
il a probablement regu la bonne commande, sinon, il doit assumer que la 
commande n’est pas la bonne et la retourner au fabricant. 

Le menuisier ne le savait pas, rnais ce qu’il a effectue est un simple ’’test 
d’hypothese”. 

4.2 Hypotheses HO et HI 

Lorsqu’ on a echantillion statistiques quelquonqe (lcs mesures des plaques), 
et on veut voir si cet echantillion verihe une hypothese donne (les plaques ont 
une taille moyenne de 50cm), on appclle ceci un test d’hypothese, on appelera 
1’hypothese qu’on veut verifier: L’hypothese nulle, note par H 0 , alors dans 
l’exemple precedent; on ecrit: 

H 0 — {les plaques ont une taille moyenne de 50cm} 

Et on notera 1’hypotese alternative, par H \, dans l’exemple precedent: 
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Hi = {les plaques n’ont pas une taille moyenne de 50cm} 

On generalise tous ce qn’on a dit dans la definition suivante: 

Definition 4.2.1. En statistiques, un test d’hypothese est nne demarche con- 
sistant a rejeter on a accepter une hypothese statistique, appelee hypothese 
nulle, en fonction d’un echantillon donne. 

Dans notre cas, on a fait une simulation d’une variable aleatoire suiv¬ 
ant une loi donnee, ce qui nous perrnet d’obtenir un echatillion de valeurs 
aleatoire issus de cette loi,dans un premier moment, on n’est pas sur si notre 
simulation est bonne, autrement dit, on n’est pas sur si L echantillion obtenu 
peut etre vraiment issu par un variable aleatoire suivant la loi en ques¬ 
tion, c’est ici oil les tests d’hypotheses entrent en jeux, pour verifier que 
L echantillion suit bien la loi en question, on effectue ce un test d’hypotheses 
d’adjustement: 

Definition 4.2.2. Soit (X\,X n ) un echantillon de variables aleatoires 
independantes qui suivent la meme loi (la loi est inconnu), un test d’hypotheses 
d’adjustement est un test qui nous perrnet de verifier si cet echantillion suit 
une loi donne. 

Remarque. II y’a d’autres types de tests d’hypotheses, par exemple le test 
d’independance qui perrnet de verifier si deux echantillions sont issus de deux 
lois independantes. 

Dans notre cas, on doit effectuer des tests d’hypotheses d’ajustements 
pour verifier si un echantillion suit une loi donne, alors l’hypothese nulle est: 

Hq = { L’echantillion suit la loi donnee} 

Et l’hypothese alternative: 

Ho = { L’echantillion ne suit pas la loi donnee} 

4.3 Test de Khi-deux (T 2 ) d’ajustment 

Definition 4.3.1. Soit X±,...X n n variables aleatoires independantes suivant 
des lois normales centrees reduites, on dit que X suit une loi de Khi-deux et 
on note : 


X ~ Si et seulement si: A" = Y^i =i X? 
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Table de valeurs des quantiles 

On aura besoin de ce tableaux pour lire des valeurs de la loi de khi-deux pour 
des degrees de libertes differentes: 


Degres de liberte Valeur du \ 2 

1 

0.004 

0.02 

0.06 

0.15 

0.46 

1.07 

1.64 

2.71 

3.84 

6.63 

10.83 

2 

0.10 

0.21 

0.45 

0.71 

1.39 

2.41 

3.22 

4.61 

5.99 

9.21 

13.82 

3 

0.35 

0.58 

1.01 

1.42 

2.37 

3.66 

4.64 

6.25 

7.81 

11.34 

16.27 

4 

0.71 

1.06 

1.65 

2.20 

3.36 

4.88 

5.99 

7.78 

9.49 

13.28 

18.47 

5 

1.14 

1.61 

2.34 

3.00 

4.35 

6.06 

7.29 

9.24 

11.07 

15.09 

20.52 

6 

1.63 

2.20 

3.07 

3.83 

5.35 

7.23 

8.56 

10 64 

12.59 

1681 

2246 

7 

2.17 

2.83 

3.82 

4.67 

6.35 

8.38 

9.80 

12.02 

14.07 

18.48 

24.32 

8 

2.73 

3.49 

4.59 

5.53 

7.34 

9.52 

11.03 

13 36 

15.51 

20.09 

26.12 

9 

3.32 

4.17 

5.38 

6.39 

8.34 

10.66 

12.24 

14.68 

16.92 

21.67 

27.88 

10 

3.94 

4.87 

6.18 

7.27 

9.34 

11.78 

13.44 

15.99 

18 31 

2321 

29.59 

11 

4.57 

5.58 

6.99 

8.15 

10.3 

12.9 

14.6 

17.3 

19.7 

24.7 

31.3 

12 

5.23 

6.30 

7.81 

9.03 

11.3 

14.0 

158 

185 

21.0 

262 

32.9 

13 

5.89 

7.04 

8.63 

9.93 

12.3 

15.1 

17.0 

19.8 

22.4 

27.7 

34.5 

14 

6.57 

7.79 

9.47 

10.8 

13.3 

16.2 

18.2 

21.1 

23.7 

29.1 

36.1 

15 

7.26 

8.55 

10.3 

11.7 

14.3 

17.3 

19.3 

22.3 

25.0 

30.6 

37.7 

a 

0.95 

0.9 

0.8 

0.7 

0.5 

0.3 

0.2 

0.1 

0.05 

0.01 

0.001 


La lecture du tableaux se fait de la maniere suivante: on choisit une valeur 
de a (a est appellc un quantile), par exemple on choisit a = 0.05, et puis on 
choisit le degree de liberte, par exemple 4, la valeur correspendante est 9.49, 
alors la lecture se fait comme se suit: 

P(X% > 9.49) = 0.05 
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Methode du test de Khi-deux 


Supposons qu’on a un echantillon de taille n issu d’une loi X qu’on ne connait 
pas, on divise cet echantillion sur des classes distinctes, la determination des 
classes depend de la loi, si c’est nne loi discete, chaqne classe va contenir un 
scule elment de l’echantillon, et lc nombre d’occurence de cet element dans 
1’echantillion va representer 1’effectif observe de cette classe, si c’est une loi 
continue, chaque classe va etre represente par un intervalle, et lc nombre 
d’elements qui apparitnet a cette intervalle var representer l’effectif observe 
de cette classe, puis on calcule les effectifs theoriques rip,, pour chaque classe 
K, , avec: 

Pi = P{X £ Kt) i £ N 

posons k = le nombre des classes obtenu”, maintenant pour effectuer lc test 
de Khi-deux, on a besoin du theoreme suivant: 


Theoreme 4.3.1. En utilisant les donnees precedentes, on a: 

C _ (nj—npi) 2 

° n ~ l^i =1 n Pi 

tend, quand n tend vers I’infini vers une loi: 


Autrement ecrit: 


hm r , 


E k 
i= 1 


(' Uj-npi ) 2 


npi 




Alors pour une valeur grande de n (Par exemple n = 50) on sait que la 
somme S n est tres proche d’une loi lc test de Khi-deux se repose sur 

ce resultat, voici la methode suivi: 


• On choisit une valeur de n (plus que la valeur est grande, plus que les 
resultats sont bien precis), dans le reste de ce chapitre, on va travailler 
avec n — 50. 


• On fait la simulation de la variable aleatoire eb question n fois, et on 
ordonne 1’echantillion obtenu d’une maniere croissante. 


• On clivise 1’echantillion obtenu sur des classes clisctinctes suivant la 
methode qu’on va expliquer, et on note par k lc nombre des classes 
obtenus 
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• Apres avoir calcule les effectifs observes et les effectifs theoriques, on 
calcule la somrne S n 

• On choisit une valeur de a, plus que cette valeur est petite, plus que le 
test est efficace (on va expliquer ceci dans un moment), pour lc reste 
de ce chapire, on travaille avec a = 0.05. 

• On lit dans le tableaux des quantile la valeur X‘l_ x pour a = 0.05. 

• Si S n ^ Xk-i 011 re j e t e l’hypothese H 0 avec une probability d’erreur 
d’au plus a , et done Pechantillion n’est pas issu de la loi en question. 

• Sinon, on garde l’hypothese H 0 , et alors l’echantillion est issu de la loi 
en question. 

Justification du procede: 

En utilisant le resultat du theoreme, on considere S'so une valeur issu de la 
loi et done: 

P(S 100 > Xl_ x ) = a 

Et done pour une valeur de a tres petite (proche de 0), Pevenement 

{^50 > 

devient presque impossible, alors si on trouve que cet evenement est realise 
(i.e S 50 > A’fc_ 1 ), cela veut dire que S 50 ne suit pas une loi A|_ 1; et done 
Pechantillion ne parvient pas de la loi en question. 

Le choix de a revient a nous, plus que a est petit, plus que Pevenement 
S 5 0 > %k-i est presque impossible, e’est pour cela qu’on choisit une valeur 
tres petite de a, cette valeur represente notre tolerance de rejeter Phypothese 
H 0 par erreur, et e’est pour cela qu’on Pappelle aussi la tolerance. 
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Part II 


Partie pratique 



Chapter 5 

Simulation des variables 
aleatoires : Partie pratique 

5.1 Inctrodution a Scilab 

Scilab (contraction de Scientific Laboratory en anglais) est un logicicl libre 
multi-plate-forme fournissant un environnement de calcul pour des applica¬ 
tions scientifiques, il possede un langage de programmation de meme nom, 
develope par 1’INRA en 1982, Scilab et souvent utilise pour des applications 
mathematiques comme il possede la plupart des entites mathematiques (ma¬ 
trices, vecteurs..etc) et leurs operations. Dans notre cas, nous allons avoir 
besoin de Scilab pour ecrire et executer les algorithmes de simulation des 
variables aleatoires qu’on a vu , done il faut lc telecharger a partir du site 
web officiel: 


www.scilab.org 

Apres avoir installe scilab, vous devez telecharger le code source du pro¬ 
gramme de simulation qu’on a develope (vous trouverz le lien a la fin de ce 
chapitre). Pour executer lc programme, il suffit de l’ouvrir a par scilab,cliquer 
‘executer’, et puis ‘sauvegarder et executer’. 
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5.2 Simulation des variables aleatoires par Scilab 

5.2.1 Loi uniforme 


On definit une fonction U qui prend deux paramters a et b. 
pour simuler la loi uniforme sur un intervalle [a, b\, cette fonction fait appel a 
la fonction rand qui suit une loi uniforme sur [0,1] et puis retroune la valeur 
centree dans l’intervalle [a,b]. 

1 function u=U(a,b) 

2 u = a + (b—a) *rand( 1 ); 

3 endfunction 


5.2.2 Loi de Bernoulli 

On definit une fonction B qui prend un parametre p, avant faire la simulation 
notre programme doit verifier que ce parametre est bien compris entre 0 et 
1 . 

Pour simuler la loi binomiale cette fonction fait appel a la fonction de loi 
unifrome sur [0,1] et puis retourne la valeur de la fonction indicatrice sur 

[°.P] 

1 function b=B(p) 

2 u = U(0 ,1) ; 

3 if u <= p then 

4 b = 1; 

5 else 

6 b = 0; 

7 end 

8 endfunction 


5.2.3 Loi de Binomiale 

On definit une fonction Binomiale qui prend deux parametres n et p. (la 
valeur de p doit etre validee par notre programme). Cette fonction alors doit 
retrouner la valeur du sommation de n lois de Bernoulli de merne parametre 
V- 

1 function x=Binomiale (n , p) 

2 x = 0; 

3 for i = 1 : n 

4 x = x + B(p); 
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5 end 

6 endfunction 


5.2.4 Loi exponentielle 

On definit une fonction expo qui prend un parametre A, et puis on applique 
le methode de fonction inverse montree dans le chapitre precedent. 

1 function x=expo (lambda) 

2 u = U(0 ,1) ; 

3 x = —log (u)/lambda ; 

4 endfunction 


5.2.5 Loi de Poisson 

Pour la loi de poisson on a etudie deux methodes de simulation, alors on a 
debni deux fonctions Poisson et Poisson2. 

Ici on va utiliser la deuxieme methode puisqu’elle est plus pratique et de¬ 
man de moins de temps, mais vous trouverz Fimplementation de la premiere 
dans le fichicr, vous pouvez essayer les deux et voir la difference entre les 
temps d’executions. 

Pour tester la premiere methode, il suffit de faire appcl a la fonction Poisson 
avec un parametre dans Scilab apres avoir execute le hchier. 

La fonction Poisson2 prend un parametre A et retoune une valeur du nombre 
de fois necessaire pour que la loi exponoticlle du meme parametre A depasse 
la valeur 1. 

1 function y=poisson2 (lambda) 

2 s = 0; 

3 i = 0; 

4 while (s <= 1 ) 

5 i = i +1 

6 s = s + expo (lambda) 

7 end 

s y = i -l; 

9 endfunction 


5.2.6 Loi normale 

On definit une fonction normale qui prend deux parametres m et u et puis 
on applique la methode de Box-Muller montree dans le chapitre precedent. 
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1 

2 

3 

4 

5 

6 

1 

2 

3 

4 

5 

6 

7 

8 

9 


function x=normale (m, u) 

ul = U(0,1) ; 
r = sqrt(—2*log(ul)); 
u2 = U(0,2*%pi) ; 
n = r*cos(u2); 
x = m + u*n; 
endfunction 


5.2.7 Loi a densite de support compact 

On definit une fonction rejetc qui prend quatre parametres: 

• a et b : les bornes du support compact de la fonction. 

• m : la borne superieure de la fonction. 

• // :une chaine de charachteres qui represente l’expression de la fonction 
en fonction de x, par exemple: " x + 1”. 

Et puis on applique la methode du rejet montree dans le chapitre precedent. 

function x=r ej et c (a , b ,m, f f ) 
deff ( ” y=f (x) ” ,”y=”+ff ) ; 
x = U (a , b j 
z = U(0 ,m) 
while (x > f ( z )) 
x = U(a,b) 
z = U(0 ,m) 

end 

endfunction 


5.2.8 Programme finale 

Apres avoir defini tous les fonctions qu’on a besoin pour travailler, void le 
plan de notre programme finale: 

• Au debut le programme doit afficher une petite documentation qui 
montre a l’utilisateur comment utiliser le programme. 

• Apres l’execution le programme demander a l’utilisateur de choisir une 

loi. 


32 


• Apres faire un choix, le programme doit demander a l’utilisateur d’entrer 
les parametres de la loi choisi. 

• A chaque entree, le programme doit verifier si les parametres entres 
sont valide, sinon il doit demander a l’utilisateur d’enter un parametre 
valide une autre fois. 

• Apres que tous les parametres sont verifies, le programme doit de¬ 
mander a l’utilisateur d’entrer le nombre de simulations (n) qu’il veut 
effectuer. 

• Puis, le programme doit effectuer les simulation en appelant la fonction 
associee a la loi donnee n fois et en affichant la valeur retourne apres 
chaque appel. 

• Apres que la simulation est fini, le programme doit demander si l’utilisateur 
veut simuler une autre loi, sinon l’utilisateur peut sortir du programme 
en entrant un charachtere specifie dans la documentation. 

Comme le code source du programme finale est long (plus de 400 lignes), 
alors pour ne pas allonger le document, on a pose le fichier du code source 
sur internet, pour le telecharger, vcuillez acceder le lien suivant: 

https://archive.org/download/simulations_ 201806 /simulations.see 
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Chapter 6 

Application du Khi-deux 


En utilisant notre programme, On fait n = 50 simulations de chaqune des lois 
qu’on va tester, comme on a deja dit, on va choisir line tolerance de a = 0.05 
(5%), vous avez lc liberte de choisr un plus grand nombre simulations et/ou 
une plus petite tolerance si vous voulez que votre test soit plus precis. 

6.1 Loi binomiale 

On veut tester si l’echantillion d’une simulation d’une loi binomiale dans 
notre programme vraimcnt suit une loi binomiale, d’une maniere arbitraire, 
on choisit par exemple une loi binomiale £>(14, |) 

Hq = { L’echantillion suit une loi binomiale -6(14, |) } 

Apres avoir fait les simulations , et calcule les effectifs observes observes et les 
effectifs theoriques on regroupe les donnees qu’on a trouve dans le tableaux 
suivant: 


Valeurs observes 

1 

2 

3 

{4,5} 

{ 6 . 0 } 

Effectifs observes n l 

14 

10 

15 

7 

5 

Effectifs theoriques npi 

7.7 

12.5 

12.5 

12.9 

4.3 


On 5 classes dans cette division, alors k — 1=4 On calcule la somme A 5 O: 

>4 (rij-npi) 2 _ 


650 - J 2 i =1 


npi 


= 8.73 


Pour a = 0.05 on lit la valeur suivante sur le tableaux des quantiles: 
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On a alors : 


Xl = 9.49 

S 50 < 9.49 

Et done one accepte 1’hypothese H 0 avec une probabilite de 5% d’avoir tort 

Conclusion: 

D’apres lc test de Khi-deux, l’echantillion suit bien une loi binomiale, et done 
cette loi est bien simulee dans notre programme. 

6.2 Loi de Poisson 

Maintenat on veut tester la loi de poisson, d’une maniere arbitraire on choisit: 

n\) 

Apres avoir effectue le precede de Khi-deux (en preservant toujours notre 
choix de n = 50), on regroupe nos donnees dans lc tableaux suivant: 


Valeurs observes 

0 

{1,2} 

Effect its observes 
Effectifs theoriques 

36 

35.8 

13 

13.4 


On a le nombre de classes k = 2, alors k — 1 = 1, On calcule la somme 
S^o'- 

Sso = EL, (jh ^- = 0.013 

Pour a = 0.05 on lit 1a. valeur suivante sur lc tableaux des quantiles: 

X] = 3.84 

Et alors : 

S50 ^ A’J 

Et done on accepte 1’hypothese H 0 avec une probabilite d’erreur de 5%. 
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Conclusion: 


D’apres lc test de Khi-deux, notre programme sirnule bien la loi de Poisson. 

6.3 Loi exponentielle 

Maintenat on vent tester la loi exponentielle , d’nne maniere arbitraire on 
choisit: 


£(§) 

La loi exponentielle est une loi continue avec une fonction de repatition F 
explicitement defini. alors pour la i’eme classe K % , pour calculer la probabilte 
Pi = P(X G Ki) Kj = [di, bi[ di, bi G M, On utilise la propriete suivante: 

P(X G = P(a < X <b) = F{b ) - F(a) 

On pent ausssi utiliser le tableaux statistique de la loi exponentielle. 

Apres avoir effectne lc precede de Khi-deux (en preservant n = 50), on 
regroupe nos donnees dans lc tableaux suivant: 


Valeurs observes 

[0,1[ 

[1,2[ 

P.3[ 

[3,5[ 

[4,H[ 

Effect its observes rq 

18 

11 

11 

5 

5 

Effectifs theoriques npi 

19.7 

11.94 

7.4 

7.1 

3.56 


On a le nombre de classes k — 5, alors k — 1 = 4, on calcule la somme 

* 5*50 • 


C _ (ni-npi ) 2 

50 A/j= 1 n p i 


3.17 


Ponr a = 0.05 on lit la valenr suivante snr lc tableaux des quantiles: 


X'l = 9.49 


Et done : 


S5 0 ^ x'i 

Alors on accepte l’hypothese H 0 avec une probabilite d’erreur de 5%. 
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Conclusion: 

D’apres lc test de Khi-deux, la loi exponontielle est bien simulee dans notre 
programme. 


6.4 Loi normale 

Pour simuler line loi normale quelquone, notre algorithmic repose sur la simu¬ 
lation d’une loi normale centree reduite, et puis effuctuer une transformation 
vers la loi normale specifie, done si la loi normale centree reduite est bien 
simulee, alors tous les autres lois normales sont aussi bien simules, si pour 
ceci qu’on choisit de tester la loi normale centree reduite : 

x ~ jV(o, i) 

La loi normale est une loi continue, alors pour diviser notre echantillion, 
on commence par la plus petite valeur observee, et puis on divise notre 
echantillion selon des intervalles semi-ouverts distincts, l’effectif observe de 
chaque classe va etre lc nombre d’elements qui appartient a l’intervalle qui 
represente cette classe, et pour la i’eme classe K tl on calcule la valeur pi = 
P(X e Kj) en utilisant le tableaux de la loi normale centree reduite. 

Apres avoir effectue le precede de Khi-deux (en preservant toujours notre 
choix de n = 50), on regroupe nos donnees dans le tableaux suivant: 


Valeurs observes 

[-2,-1[ 

[-i,o[ 

[0.1[ 

[1,2[ 

Effectifs observes n t 

11 

16 

17 

6 

Effectifs theoriques npi 

6.8 

17.1 

17.1 

6.8 


On a le nombre de classes k — 4, done k — 1 = 3, on calcule la somme 


S^o- 


_ V^4 (m-npi) 2 _ 


£50 - Di=1 


npi 


= 2.75 


Pour a = 0.05 on lit la valeur suivante sur le tableaux des quantiles: 


Et alors 


A 3 2 = 7.81 


57 n ^ X4 


150 ‘' l 3 

Et done on accepte Fhypothese H 0 avec une probabilite d’erreur de 5%. 
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Conclusion: 

D’apres le test de Khi-deux, notre programme simule bien la loi normale 
centree reduite, et par consequence toutes les autres lois normles . 
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